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Abstract 

Using large-scale real-time lattice simulations, we calculate the baryon asymmetry generated at 
a fast, cold electroweak symmetry breaking transition. CP-violation is provided by the leading 
effective bosonic term resulting from integrating out the fermions in the Minimal Standard Model 
at zero temperature, and performing a covariant gradient expansion [l|. This is an extension of the 
work presented in [2J. The numerical implementation is described in detail, and we address issues 
specifically related to using this CP-violating term in the context of Cold Electroweak Baryoge- 
nesis. The results support the conclusion of [2], that Standard Model CP-violation may be able 
to reproduce the observed baryon asymmetry in the Universe in the context of Cold Electroweak 
Baryogenesis. 
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I. INTRODUCTION 



At its inception, Electroweak Baryogenesis was an attempt to explain the observed baryon 
asymmetry of the Universe by processes originating from Standard Model physics only |3|, |4|. 
Although essential ingredients, C-, P-, CP- and baryon number violation are present in the 
Standard Model, the detailed quantitative implementation of the scenario has encountered 
two stumbling blocks: At experimentally allowed Higgs masses, the electroweak transition 
is an equilibrium cross-over [5], rather than the required out-of-equilibrium phase transi- 
tion; and at electroweak-scale temperatures, the effective CP-violation is much too small to 



produce the observed asymmetry 



The way to alleviate these problems has traditionally been to embed the Standard Model 
in a larger theory, either by exten ding th e Higgs sector (see, for instance |9|), or by imposing 



full-fledged supersymmetry (see [10Ml5f| for recent developments). This allows for a first 
order electroweak phase transition, while introducing additional coupling constants in the 
Higgs sector, which if taken complex may give rise to in principle an arbitrary amount of 
CP-violation. 

Standard Model CP-violation originates from the CKM-mixing in the fermion mass ma- 
trix, but is often phrased in terms of effective bosonic terms in the action, appearing as a 
result of integrating out the fermions in the path integral. The strength of CP-violation is 
encoded in the coefficient of the CP-violating terms, and is proportional to 

J (m 2 t - ml) (m 2 t - m 2 u ) (m 2 c - m 2 u ) (m 2 b - m 2 ) (m 2 b - m 2 d ) (m 2 s - m 2 d ) (1) 

in terms of the quark masses m^ and the Jarlskog invariant flf] ]. 

J~3.1xl(T 5 . (2) 



At high temperatures T > m i} the coefficient is known to be very small oc J(mj/T) 12 and 
electroweak baryogenesis at these temperatures require additional sources of CP-violation. 



At zero temperature, the situation changes radically. In 17H22| two separate methods 
have been developed to calculate the CP-violating part of the effective bosonic action at 
zero temperature. Integrating out the fermions, we have 

4>,A] = J V(j)VAV^>V^e lS ^^^ A] = J V( j)VAe iS[<t '' A]+iTrLosD[<t '' A] (3) 



2 



with D the fermion propagator and <f> and A the Higgs and gauge fields, respectively. In 



both approaches, the fermion contribution is then expanded in covariant gradients 



leading order (order 4 in derivatives), CP-violation is absent [21], but at next-to-leading 



4J]. At 
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order (order 6), CP- violating terms have been reported in the P-odd [1] and P-even 
sectors. The authors of [22| also finds vanishing CP-violation at order 6 in the P-odd sector 
in conflict with [1]. 

The crucial point is that at zero temperature the coefficients of these terms are only 
suppressed by J, but not by powers of rrii/v, v the Higgs vev, as one may have expected. 
This opens up the possibility that Standard Model CP-violation may be responsible for the 
baryon asymmetry, if baryogenesis took place at low enough temperature. 

Cold Electroweak Baryogenesis is realized if the temperature is well below the electroweak 
scale when the electroweak transition happens 23M26l| (further developed in 27H3l|). In the 
Standard Big Bang scenario, inflation and reheating occur at very high temperature, and in 
the Standard Model the electroweak cross-over takes place around T = 100 GeV. However, 
if Higgs symmetry breaking is triggered not by the change in the finite-temperature poten- 
tial but by the coupling to another field, the transition can be delayed until cosmological 
expansion has cooled the Universe far below the electroweak scale. Different mechanisms 
to produce this transition have been proposed including low-scale hybrid inflation 25 
two-step symmetry breaking and a specific variant of a first order phase transition 

In terms of the Higgs potential, we write 

^(0)=/i e 2 ff (t)0V + A(0V)2 ) (4) 

with the electroweak transition happening when /i 2 ff (£) = 0. For the purpose of the present 
work, all we need is that the Higgs mass parameter flips sign quickly, so that to a good 
approximation 

fi 2 eS (t < 0) = /i 2 , fi 2 cS (t > 0) = -/i 2 , (5) 

i.e. an instantaneous quench. The system goes through a spinodal (or tachyonic) instabil- 
ity, where field modes with k < /1 grow exponentially until Higgs self-interactions become 



28], 
33). 



significant 



34j. This stage of (p)reheating leads to large field occupation numbers in the IR 



35M40j . making the dynamics effectively classical 



36 



38|. 



In the following, we will assume that Standard Model CP-violation manifests itself at 
leading order as the CP-odd term of lj], with some coefficient k cp , and that the asymmetry 



scales linearly in k cp (this is certainly true for small values of k cp ). In numerical simulations 
of a cold electroweak transition, we will calculate the baryon asymmetry generated for a given 
k cp . This will allow us to estimate a value of k cp , above which Standard Model CP- violation 
can accommodate the observed asymmetry. Awaiting the resolution of the discrepancy in 
Q, Q]> we will consider the specific CP- violating term used here as representative of all 
P-even and P-odd terms at NLO in the covariant gradient expansion. The important and 
common feature to both results is that the coefficients of these terms are not suppressed by 
powers of the Yukawa couplings (rrii/v). 

We will follow the dynamics of the system starting from a vacuum (T = 0) initial state, 
and calculate the baryon number asymmetry produced in the transition under the influence 



31 



411 ]. except 



of the CP- violating term. The implementation closely mirrors the work in |30l . 
that the CP- violating term is different (and much more complicated). 

Preliminary results were presented in where we found that Standard Model CP- 
violation as represented by the leading bosonic term can indeed account for the observed 
baryon asymmetry. The present paper is an update confirming this conclusion with signifi- 
cantly more statistics and providing the numerical details of the simulations. 



II. SU(2)-HIGGS MODEL WITH CP- VIOLATION 



Serving as an approximation to the full Standard Model, we consider the SU(2)-Higgs 
model with CP-violation, with continuum action 



d 4 x 



1 



^Tri^F^ + (Drf? Dy + /i e 2 ff (t)0V + A (0V) + V - S CPfi 



(6) 



(f)(x) 



(7) 



in terms of a complex scalar doublet and an SU(2) gauge field 

01 (x) 
4> 2 {x) 

where a a are the Pauli matrices. We have F^ v = d fM A u — d v A^ — i[A^A v ], = 
(dfj, — iAfj) <fi. The vacuum particle masses are m w = gv/2, m H = a/2/x = \/2Xv in terms of 
the Higgs vacuum expectation value v = fi/y/X, with Vq = /i 4 /(4A) normalising the energy 
to zero in the vacuum. 
We also define the useful 



and note that 

= (e±f = e* (e°) f = (e°) T = e°. (9) 
We want to write everything in terms of the Higgs matrices, with {45] 

$ = ( <M^) 0i(^) ] = I ^0*0 j ^ 10 ^ 

\^ 2 (z) 2 (x) J \-<j>*(x) fa{x) J ' 

or 

4>(x) = ia 2 0*(x). (11) 



We have 



i 2 | = -Tt^ + ^, (12) 



and for some matrix B 



•*(«)*(«, „)•<») = I ^ \. (13) 

^<ft(x)B(x,y)<j)(y) <p ] (x)B(x,y)(j)(y) 

We then define 

B° = Tre°^\x)B(x,y)^(y) = ^(x)B(x,y)<l>(y) - 4>\x)B(x,y)4>(y), 

(14) 

5+ = Tre + &(x)B(x,y)$(y) = 4t(x)B(x,y)ij>(y), (15) 
B- = Tre~& (x)B(x,y)<l>(y) = ft (x)B(x,y)(j>(y). (16) 

We are particularly interested in the case when B is composed of covariant derivatives D, 
and we define the following: 

B -> Dj.foy) : C£'° = Tre ± '°^(y)D fl (x,yMx), (17) 
B -> D v (x,y)D x (x,y) : L>± = 'B: e ± $t(y)D I/ (ar,j/)D A ( a r,j/)$(ar). (18) 

We note the important relations under complex conjugation 

(Cj)* = -Cy, (Cj)*=-C;, (e^D*)* = -e^Zft, (19) 
which follow from the definitions of the covariant derivative, 

e ^Xa F u\ = e ^\a^ D ^ D ^ = 2ie^ a D v D x , (20) 



The CP-violating term can be written [l| as 



in 



S, 



CP,6 



d x ■ 



IK 



£ ,uxa i c o c+c-c- + c;c+c~ D uX 



K C aC+C+ + CtC- CI )D uXj . 



(21) 

where <p\ will be defined below. Note that the two terms are each other's complex conjugates, 
but with a minus sign [46]. Using v = 246 GeV, raw = 80.4 GeV and m c = 1.3 GeV, 
J = 3.1 x 10~ 5 , we have 



re ~ 1.3 x 1(T 4 k cp . 



(22) 



The Standard Model, zero temperature, value for k cp was found in [l| to be 9.87. 



A. The cut-off A 

The CP-violating term is divergent at points in space where <f>*^> = 0, where the gradient 
expansion breaks down. Although rare once the electroweak transition has begun, at any 
given time during the simulation, there will exist small regions of space where the CP- 
violating force is "unphysically" large. Apart from discarding the approach altogether, the 
solution is to identify these regions and cut them out of the CP-violating dynamics. This 
can for instance be done by imposing a threshold A, and set the CP-violating force to zero 



whenever ^J^<j){x) < A. Rather than a step function, in practice we will we introduce the 
cut-off Higgs field (j>\ 



"> U = TTUTTT^i ( 23 ) 



J A 



c(0t0 + A 2 ) 



fixing the number c by 



l + ^ : 

V 



(24) 



so that 



0V = t^ cw 1 + a 2 ) = ^2- (25) 

We emphasize that the cut-off A is in no way related to a momentum cut-off. It is simply 
a threshold which moderates the divergence in the CP-violating force near small |</>(x)|. We 



also note that the introduction of such a threshold means that the result of the simulations 
will be a lower bound on the generated asymmetry from the Standard Model, and there will 
be a dependence of the result on A. In 2J we used A = 50 GeV. Below, we will study in 
detail the dependence on A. We note that A G [0, v/y/2] = [0, 174] GeV. 

III. LATTICE IMPLEMENTATION 

We here explicitly write down the lattice action and equations of motion for including 



30 



31 



41], and may be skipped by readers who 



CP-violation. This runs along the lines of 
are not interested in the numerical implementation. 

On the lattice, we perform a rescaling of the Higgs field 

$ V\a x $, (26) 

We use the usual lattice derivatives 

apD^ = U Xil & x+ll - $ x , a„D'^ x = <& x - U^n&x-p. (27) 
The symmetrized derivatives are 

d; = 1 (Dp + , d% = d;d° u . (28) 

Then let us define a slightly different lattice version of our factors for \x 7^ 

4>\{x)C;{x) = 2a^Tr[e%D s ^ x ] = Tr[e%+ (lJ^ x+ , - C^ M $_ M ) ] , i = +, - 0, 

(29) 

and a special version for the timelike factor /i = 0, 



(30) 



[x)C*(x) = Tr[e* (*1 A>*. + ^ x+0 D^ x+0 ] 

= Trie* ($i(U , x <S> x+0 - 9 X ) + $t +0 ($ x+0 - 1 = +, -, 0. 



This is necessary to avoid the equations of motion from becoming implicit two steps ahead 
in time. In this way they are implicit one step only. We note, that 

Tr[e 4 $t$,] =0, i = +,-0. (31) 



We also define, for //, v ^ 

Again, special rules apply for the timelike derivatives 



(32) 



= Tr 



+ (^l+o(^,*+0^x +V +0 " C/ t ,x^,^x+, - C/^+q-^x+O-, + ^.x^x-^x-,)) ] , 



(33) 



and 



0i(a;)^ o = Tr e ± ($l(U UtX U , x+u $ x+0+u - U n , x $ x+V - Ut tX _ v U , x - v $ x - v+0 + Ut tX _ v $ x - v )) 
(®l+o(Uv, x +o®x+v+o - U UjX+0 Ul x+u $ x+u - ul >x _ v+0 $ x - v+0 + uI >x _ v+0 Uq >x _ v ^ x - u ^ 



+ 



We note that compared to our continuum notation we have the substitution rules 



(34) 



-ii \ continuum 



<f> 2 A (x)(CX mce 2a M (C;) c 
<Pi(x)(D l uX ) latUce <-> Aa v a x {Dl x ) contmuum . 



(35) 
(36) 



Then the lattice contribution to the action is 



SCP,6 - ^2 ^ UXa ~jJ~ \ C l,x D t\,x ( C a,x C a,x C a,x + C o-,x C a,x C a,x) ~f 
x,t -^ V a ° 

1/7O n- (n- n+ n+ 1 /^+ ^- ^+ \ ° 

^ u,x v\,x \ <r,x a,x a,x ~ cr,x a,x a,x) 



(37) 



with the dimensionless 



3JK CP (alm%/4) 



2 ie n 2 a 2 x m 2 c 



(38) 



In the following we will use the shorthand 



cfivXcr — r uv\cr u x 



(39) 
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The complete lattice action of the SU(2)-Higgs model with CP-violation then reads, 

x,t n ^ ' 

+/^Tr[(f/ ,A + o - $(z)) f (U ^ x+0 - $(x))] 

E ~Tr[(U niX $ x+n - $(x)) f - $(x))] 

n 

-fo QMt$ B - vl t ) ] - S CP>6 , (40) 

where by matching to the continuum theory, we have in addition to (|39|) 

at _ 4 a * os _ 4 a t R t _ la x os _ la t o _ la t 2 _ K^) 2 
Pg ~ ~~o ' PG — ~~o ' PH — T ' Ph — T 5 PR _ T 5 u Ut — Z • 

g at 9 a x a at Xa x \a x 4 

(41) 



A. Higgs equation of motion 

The Higgs equation of motion, in the A = gauge, reads 



1 SS, 



CP,6 



(42) 
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with 




Z Xn P« c ( 6 



x 2Im(c (y))$. 



-2^° [($, +0 - $ ,K(y) + - <VoK0/ " 0)] 
+2<T° [(<Vo - S^G/) + ($„ - $ H )Bf 2 (y - 0)] 
+2^° [(^ +0 - $ w )SS(y) + (*„ - ^-0)^2(2/ " 0)] 
-2(1 - O - t^^v-/.) 5 3 (y) 

+2(1 - O - C/^-^,-.) flf 2 (y) 

+2(1 - O (Ua,^ - 5f 2 (y) 
+2(l-0(l"0^a(») 

X (^.x^A.x+^x+^+A + Ul,x-vU\,x-v-\®x-"->> ~ Ul,x-vU\,x-v®x-v+\ ~ U VjX U{ x+v _ x <b x+v -\ 

+2<T°(1 - <5 A0 ) x (f/^+o^+o+A - U Ky % +x - U[ y _ x+0 ^ x+0 + U{ y _ x ^ x )B^{y) 
+2<f"°(l - <5 A0 ) x {U x , y % +X - t^ >tf _oVo+A - + <y-A-o*v-A-o)*£(v - 0) 

+2(1 - cT> A0 x (C/a,^, + o + , - U Uty % +u - Ul y _ u %_ u+Q + ^.^V^i* (v) 
+2(1 - 5 M )5 X0 x ([/,,,$,+, - U v , v * v -*+„ - Ut >y _ u %-„ + Ul^ u %^ Q )B^{y - 0)). 



(43) 



We have defined 



cM = K(c:c-c- + c-c^c-)Dt x ) 



(44) 
(45) 
(46) 
(47) 
(48) 



c„(y) = (C*C-C- + C;C+C~)Dt x ) 

c.(y) = (cJC7-(C7-£)+ +q+i); A )) w , 




c Q (y) = (C;[2C CT -C Q + L»; A + C:C~D- X + C"C7-£>+]) 

C A (y) = K(C:C-C- + C-CjC-)) y , 
„.u _ (^ 1 Re(c a (y)) - i(7 2 Im(c a (y))) 



(49) 



(50) 
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B. Gauge field equation of motion 



For the gauge fields, we use the rules 

j^-U^ = iSyU^l, j^Ul x = -iU u , y S2a b 5;i (51) 

The matrix S£ is unknown, but we will not need it, as it cancels out at the end of the 
calculation. We assume it is invertible. We go to temporal gauge A = at the end. We 
introduce the "electric" field asj^ 

El x = --Tr^X.+ol- (52) 

We then have 

^G^0^n,y ^~ D m Tr [i(T Um^yUn^y+mUm^+nU^y] 

m 

-/3 s H Tr[(U n>y % +n - *J ia a $ y ) - (S b a )- l6 ~^ = 0, (53) 

0/i y,n 

where we have introduced the adjoint (backwards) covariant derivative 

(D'm)xyFy = F a x - \^[U x . m ^Ul_ m ^ b ]F b x _ m . (54) 

We need to define 

F b = ^ y ia b U n ^ y+n , (55) 

Hl b v = $>l_ i<T b U n , y $y +n , (56) 

H2 b y = ^ y+0 ia b U n , y % +n , (57) 

HSl = ^ y+n _ Ul y ia b %, (58) 

H = ^ y+n+0 Ul y ia b %, (59) 

Gl b y = & y ia b U n , y U Ky+n $ y+n+x , (60) 

G2 b y = ^ y+n Ulja b Ul y _ x %_ x , (61) 

G3l = $l +n Ul y ia b U x ,y% +x , (62) 

G4 b y = $X[/„/i J+n _ A V-A, (63) 
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D„\ (C^C a C a + C a C^C a ) 



(64) 
(65) 
(66) 
(67) 



BUv) 



cl(ctc~c- + c;c+c-) y , 

{ia l Re{c a {y)) - ia 2 lm(c a (y))) 



(68) 



B$(y) 



(ta 3 Re(c a (y))) 



(69) 



The last term in fl53l) reads, 




-8^[Tr(B%(y)-B$(y + n))F;] 

+5 na l?r{BUy)-BUy + n))Fl] 

~C [Tr«(y) - B#(y + n))(^2{ - HA b y )] 

~C [Tr(5r 2 A (y + n - 0) - 2#(y - 0)){Hl h v - #3»)] 

+<y w (l - 5 0A ) [Tri?r 2 A (y)(GlJ - G4j) + Tr^ 2 A (y + n)(G3 b y - G2 b y ) 

-5™(1 - 5 0X ) [TrB$(y + X)G3 b y + TrBg(y - A)G2 6 J 

-8™(1 - 5 0X ) [TrB$(y + n + X)Gl b y + TrB^(y + n- X)GA b y }). 



C. Gauss law 

Gauss Law is the equation of motion (or rather, constraint equation) resulting from 
variation with respect to A . It means that the quantity 



(70) 



-1 




(71) 
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is constant in time. We first define 



F b = 

y 




(72) 


m b v = 


H mb (Ux, y+ o% + o+x - Ul y+0 _ x <f> y+0 _ x ) , 

\ / 


(73) 


H2 b y = 




(74) 


m\ = 




(75) 


HA b = 

V 




(76) 


cM = 


[(c+c- + c;c+)c-D+ x ] y , 


(77) 


c*{y) = 




(78) 


c a \y) — 






Cu\(y) = 




(80) 


BUy) = 


(icr 1 Re(c a (y)) - icr 2 Im(c a (y))) 

#(v) 


(81) 


B a ,(y) = 


(^ 3 Re(c a (y))) 
</> c 2 (?/) ' 


(82) 



and then we have 

= Wy~f" Xa ( ~ 2Sfl0TrF S B " to + 2^°1VF^r 2 (y) + 2^ TVFX(y) 

<T°Tr(/nJ + H%)B$(y) - 5 X0 TrHA b y B^(y + v) + ftriJ^fe - i/) j . 

(83) 



IV. NUMERICAL PROCEDURE 



The equations of motion are implicit in time, but because of the specific choice of time- 
discretization, there are only nearest-neighbour couplings in time, while there are next-to- 
nearest neighbour couplings in space. This makes it possible to solve the equations of motion 
by iteration in the CP-violating force term. Because the size of the CP-violating term is 
small, this converges under certain conditions on k cp and dt. 

We use a lattice of V = L 3 sites and a lattice spacing a x , with 

a x ^H = 0.35, mn/mw = 2, L — 64, g 2 = 4/9, dt = a t /a x = 0.0125, 

(84) 
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and we run the simulations until the system settles, m#t stop = 30. We will choose kcp 
to get an observable signal at a given cut-off, and then scale back to a common k cp value 



assuming a linear dependence 48j . The combination of large lattices, small timestep and the 
need to iterate makes this a numerically very heavy problem, requiring O(10 6 ) CPU hours 
on state-of-the art computer clusters. 

A. Initial conditions 

We generate an ensemble of initial conditions, with vanishing gauge field = 0, and the 
Higgs field reproducing the quantum vacuum in the pre-quench potential V(4>) = +fi 2 ' 



(4)^0 = 77^, WlW^o ^^^ 2 , \k\<H (85) 

Only the unstable Higgs modes are initialised, in order to mimic that the quantum modes 
grow due to the spinodal instability, and that we only keep the modes that thus become 



classically large 35l437l|. The spinodal instability affects all the modes with /q 2 at < Xv 2 , which 



means sets of integers k x , k y , k z , satisfying 

{ -^^>kl + kl + kl (86) 

For our choice, we have ~ 80 unstable modes, enough that the IR dynamics is well repre- 
sented. 

The initialisation is generated by Monte-Carlo sampling in order that the total charge on 
the lattice is zero. However, the local charge density is non-zero, and the gauge momentum, 
the electric field E®, follows from the Gauss constraint, given the Higgs background, 

IP, 1 

%K = i^Tr[(* x+0 - S,) WJ. (87) 

Pg 

The solution to the Gauss constraint makes the approximation that k,cp — 0. Whereas in 



271 . l30l . l3l| the CP- violation vanishes when = 0, in the present case, there is a nonzero 



term which includes the electric field E, namely 



6 X0 TrH4 b y B$(y +v) + S X0 TrH3 b y B^(y - u) 
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(88) 



Solving the full Gauss constraint now becomes a non-linear minimization process, which we 
do not attempt here. We note that we do not observe any effect of this approximation. The 
asymmetry generated does not arise in the early stages of the simulation, but after the initial 
Higgs roll-off. In the absence of CP-violation, the Chern-Simons number averaged over the 
initial ensemble is indeed zero. 



B. Observables 



We will monitor the Chern-Simons number 491 ]. 



I * 

Ncs(t) - iV cs (0) = y^J2J2 [^W^Wl > (89) 

x 

the Higgs winding number, 

x,ijk 

(90) 

with M = $/|0|, and the average Higgs field, 

^ = ^J2l Ti ^- ( 91 ) 

X 

Through the anomaly equation, the Chern-Simons number (1891) is directly related to the 
baryon number of fermions living in the gauge field background. However, in the vacua the 
Chern-Simons number is integer and equal to the Higgs winding number. The latter is a 
much cleaner observable since it is always integer {50I and settles early on. We therefore use 
that at asymptotically late times 

B(t) - B(0) = 3 (N cs (t) - N cs (0)) = 3 (N w (t) - N w (0)) . (92) 

In addition, in order to minimize the statistical noise, we average over a strictly CP- 
symmetric ensemble of initial field configuration, including for each random initial con- 
figuration its C(P) conjugate. In practice, this amounts to running each configuration first 
with +/3 K , calculating the Chern-Simons N cs and the Higgs winding number iV w . And then 
running again with — (3 K , calculating —N cs and —N w . In this way we can quantify the 
CP-asymmetry configuration by configuration, by calculating the integer 

AN w = N w (P K )-N w (-f3 K ), (93) 
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and perform the statistics on these, rather than the N w individually. This has the advantage 
of reducing the statistical errors, since we note that in practice AA W is either 0, ±1 or ±2. 
A more detailed exposition of this procedure can be found in [30 1. 



V. RESULTS 



A. Size of the force 



To gather insight into the size of the CP-violating effect and the impact of the cut-off, it 
is instructive to first perform simulations using CP- conserving dynamics, while at the same 
time calculating the CP-violating force from the gauge equation of motion. Writing the 
discretized equation (!70|) as 

E%(x,t) = E°(x,t - dt) + 5E a f(x,t) + (3JE^(x,t). (94) 

in terms of a CP-symmetric SE^° and a CP-breaking force component P K 5E^' 1 we average 



the latter 



5l| 



i^T/W^M), (95) 

X 

for all values of a = 1, 2, 3, /i = 1, 2, 3. We vary the cut-off y/cA through the values 

v^A = {0, 1, 3, 10, 30, 90} GeV. (96) 



Fig. [T] shows the average force in time. The black dashed line is (for clarity, a rescaled 
and inverted version of) the Higgs field squared, averaged over the lattice (l9~Tj) . cf^cf) "rolls 
off" the potential hill from zero to some finite value, around which it oscillates. This value is 
the broken phase expectation value at finite temperature, in this normalisation about -4.9. 
For comparison, the zero temperature vev would correspond to -5. 

The oscillations of the averaged Higgs length disguises that each field configuration is 



highly inhomogeneous. As was shown in 



29 



there are many zeros initially and a few 



whenever the average Higgs amplitude is small. The distribution of zeros is in turn strongly 
dependent on the value of the Higgs mass [291 ] . 

At such nuclei, where the Higgs field length is very small, we expect an amplification of 
the force term. Indeed we see in Fig. [1] that the force is largest in peaks corresponding to 
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FIG. 1: The logarithm of the average force term in a run with dt = 0.05, n = 10, cut-offs (black), 
1, 3, 10, 30, 90 (green) GeV. Superposed, the average Higgs field; to fit into the picture, multiplied 
by -5 (black dashed). 



the extrema of Higgs oscillations, with a small delay compared to the actual minima. This 
feature corresponds to the observation that the maximum number of Higgs zeros also lags 
behind the minima of the oscillations of the averaged Higgs length 29]. 

Without a cut-off (black) the averaged force is very large, especially initially. By studying 
the distribution of the force in space, we found that the average is indeed dominated by a few 
isolated points. The force is also strongly peaked in time, following the Higgs oscillations. As 
we introduce a cut-off and increase it, as long as a/cA < 10 GeV, we see no major difference 
(red, blue lines). But for larger cut-off (10 — 30 GeV, magenta, orange lines), the peaks are 
cut down. Finally at a cut-off of 90 GeV, the peaks are almost completely gone. 

On a technical note, the light green dashed line is a representation of the value above 
which the iteration algorithm has been seen to not converge. At the values of k cp used, we 
had 4-6 iterations per timestep. 

The force is strictly proportional to the CP violation prefactor, k cp . To get maximum 
asymmetry for a given cut-off, we should use a k cp bringing the force close to the dashed 
green line. On a lattice of the size used here, a few field configurations will then produce an 
asymmetry, and we should interpolate back to the "physical" value of k cp = 9.87 to find the 
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FIG. 2: Left: The average Higgs length (black), the Chern-Simons number (red) and the Higgs 
winding number (blue) for a pair of trajectories with k cp = ±50. The cut-off is yfcA = 100 GeV. 
Right: The average Higgs length (black), the Chern-Simons number (red) and the Higgs winding 
number (blue) for the same pair of trajectories as in on the left with ac cp = ±50. The cut-off is 
now y/cA = 50 GeV 

actual baryon asymmetry. We will do so assuming that the dependence of the asymmetry 
on k cp is linear. We note that the peak values of the averaged force roughly scales as A~ 4 . 
Still, this is an averaged quantity, and it may be that at a nucleus the scaling is stronger. 

B. Single trajectories 

Let us consider the effect of CP-violation on a single pair of trajectories, an example of 
which is shown in Fig. [2] (right). The cut-off is y/cA = 100 GeV, and we show the Higgs 
length (black, dashed line), the Chern-Simons number N cs (red) and the Higgs winding 
number N w (blue). Shown are two trajectories, corresponding to /t CP = ±50. The Higgs 
field falls off the potential hill, and starts oscillating around the minimum. Meanwhile the 
Chern-Simons number and winding number of the two configurations move closely together, 
except for a small wobble in the Chern-Simons number just after the first minimum of the 
Higgs oscillation. As we have seen in the previous section, this coincides with the occurrence 
of Higgs zeros and a peak in the CP- violating force. 

In Fig. |2] (right) we show the exact same configuration, but run with cut-off y/cA = 
50 GeV. Now the discrepancy between the two sets of lines is much bigger and present in 
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Higgs Nw + 




FIG. 3: The average Higgs length (black), the Chern-Simons number (red) and the Higgs winding 
number (blue) for the same pair of trajectories as in Fig. 3 and 4, with k cp = ±800. The cut-off 
is y/ch. = 100 GeV. We have an asymmetry, AN W = +1. 



both the Chern-Simons number and the Higgs winding. In this case, the CP-violating force 
is not large enough to drive the Higgs winding to different final integers, and so AN W = 0. 

In Fig. [3j we show again the same pair of trajectories but now for k cp = 800, cut-off 
\/cA = 100 GeV. Indeed, we do now have an asymmetry, AiV w = +1. We also see that 
the Chern-Simons number is smooth, as a result of there being very little UV noise. All 
the power is in the IR tachyonic modes, until times one or two orders of magnitude longer 



than considered here 



38 



a-|40(. There is therefore no need to "cool" the gauge field in order 
to obtain a reliable value for the Chern-Simons number, in contrast to simulations of the 
(classical) equilibrium sphaleron rate (see for instance 42|). The Higgs winding number is 
expected to be discrete up to lattice discretization artefacts. This is indeed the case at late 
times, and we consider the transition over and done by time m#t = 30. Most of the action 
seems to take place at the first minimum in the Higgs oscillation (see also below) around 
rant = 10. At early times, the Higgs field is full of zeros and the Higgs winding is ill-defined, 
hence the wild behaviour up to the point when the Higgs field has a achieved a significant 
fraction of its broken phase value. Finally, the Chern-Simons number and Higgs winding 
appear to roughly move together and at the same time. The Chern-Simons number is not 
constrained to be integer, except in the vacua. Also only dynamical constraints force it to 
follow the Higgs winding, thus minimizing the covariant derivative term in the energy. This 
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FIG. 4: Left: Ensemble averages of the Higgs length (black), the Chern-Simons number (red) 
and the Higgs winding number (blue) with ft CP = ±250 and the cut-off is yfck = 75 GeV. Right: 
Ensemble averages of the Higgs length (black), the Chern-Simons number (red) and the Higgs 
winding number (blue) with k cp = ±800 and the cut-off is \fch = 100 GeV. 



also makes the Higgs winding number the cleaner observable. 



C. Ensemble averages 

We perform simulations of an ensemble of (pairs of) initial configurations, and determine 
from each pair the integer AN W , which we then average over. Fig. H] show the average Higgs 
field (full black), the Chern-Simons number (red) and the Higgs winding number (full blue) 
for a cut-off of 75 GeV and 100 GeV. The simulations were done with k cp = 250 and 800, 
respectively. N cs and iV w seem to move at the same time, and also in the ensemble averaged 
observables there is a clear asymmetry, largely created during the first minimum of the Higgs 
oscillations. This is due to the appearance of actual Higgs zeros, where winding happens 
easily, but may also be because the CP- violation is larger in the presence of such points. In 
fact, since a factor of 100/75 in y/cA requires a factor of 800/250 in k cp to get roughly the 
same asymmetry, the result is certainly sensitive to the cut-off. More about this in the next 
section. 

The dashed blue line is the average of AN W , with the dotted lines the error bars on that 
number (standard deviation). The results are presented in Table IVCl 
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TABLE I: The cut-off, k in the simulation, the number of AiV w = ±1 configurations, total 
number of configurations N tot , the corresponding asymmetry at k cp = 9.87, and the required k cp 
to match observations, all assuming a linear dependence on k cp . The last line is an extrapolation 
to y/cK = 246/\/2GeV assuming an exponential depencence on the cut-off. 



D. Cut-off dependence 

We can convert the average the Chern-Simons number to a baryon asymmetry, by writing 

M = *J&± = ™ (97) 
The density of photons is given as 

1 2vr 2 



9*T T 3 eh , (98) 



' 7.04 45 

in terms of the number of relativistic degrees of freedom g* = 86.25 and the reheating 
temperature T re h. This is in turn given by distributing the Higgs potential energy on the 
degrees of freedom, 

Vo = ^ = ^9*T* h . (99) 

Fig. [5] is the main result of this paper, and shows the dependence of the final baryon 
asymmetry as a function of the cut-off a/cA, rescaled to k cp = 9.87. There is a fairly 
strong dependence, best approximated by an exponential ub oc exp(A/cA/20 GeV), but also 
consistent with a power law % oc (-y/c A) -4 ' -5 ' -6 . 

Clearly, without knowledge of the actual effective cut-off, we cannot make a definitive 
prediction about the precise baryon asymmetry. However, for the entire allowed range of 
cut-offs (left of the vertical dashed line), the asymmetry rescaled to k cp = 9.87 is at least 
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FIG. 5: The dependence of the final baryon asymmetry n^/n 7 on the cut-off \fck. The purple 
lines are exponential forms (note the log-scale), while the red, green and blue lines are power laws. 
The horizontal pink line is the observed asymmetry, the vertical black-dashed line the maximum 
cut-off v/y/2. 



an order of magnitude larger than the observed value of 

6 x 1(T 10 , 



43] 



n B 

n. 



(100) 



7 



represented here by the pink horizontal line. An alternative way of stating the result is that 
for any cut-off, there is a k cp for which the asymmetry is equal to the observed one; and 
this k cp is always less than the zero-temperature value of 9.87. See also Table IV CI 



VI. CONCLUSION 



We have performed complete non-perturbative simulations of a cold, fast electroweak 
transition in the presence of Standard Model CP- violation, represented by f )37|) . A non- 
physical effect of Higgs zeros is dealt with by a cut-off A, and we calculate the resulting 
baryon asymmetry and its dependence on A. We again stress that this is not a momentum 
cut-off (see section III Ap . When taking the coefficient k cp to have the zero-temperature 
value reported in [l|, we find an asymmetry larger than the observed one by a factor of 
between 10 and 10 4 , depending on the cut-off. 
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The total computer time used for the results presented here was approximately 10 6 single- 
CPU hours on a state-of-the-art computer cluster. Although it may be possible to speed 
up the code somewhat, this makes large-scale parameter sweeps in Higgs mass, quench time 
and k cp prohibitively computer intensive. The lattice size was chosen so that each field 
configuration could fit on a single CPU, making parallelisation unnecessary. We believe that 
the lattice is big enough to correctly reproduce the IR dynamics (the spinodal instability), 
in which case a scaling of lattice volume is equivalent to the same scaling of the number of 
configurations [52 ] . 

There remains several avenues for improvement: Firstly, the cut-off introduced shields us 
from the breakdown of the gradient expansion near Higgs zeros. It would be interesting to 
attempt a calculation of this effective cut-off, or even better find a suitable (non-gradient) 
expansion near such nuclei. It would be interesting to see whether the same CP-violating 
term also appears at leading order in such an expansion. Secondly, so far the coefficient 
k cp has only been calculated at zero temperature, i.e. with zero external momenta running 
into the fermion loop. Although the initial condition before the electroweak transition 
is cold, during the spinodal roll-off, Higgs and gauge fields are far from equilibrium but 
also far from zero temperature. Hence, it would be appropriate to calculate k cp in the 
background of "post- spinodal" gauge-Higgs fluctuations. This is likely to be smaller than at 
zero temperature, but how much depends on the details of the transition, and in particular 
the speed of the quench. 

Ultimately, one would like to include the fermion fields themselves in the dynamics, with 
the full CKM matrix. This implies implementing fermions on the lattice in real-time; on 
the other hand there will then be no problems with Higgs zeros or numerically complicated 
bosonic terms. 

In simulations with the leading order CP- violating term 53], sufficient baryon asymmetry 
implies a bound on the quench rate \3l\ . 

i*4W (101) 

2/x 3 dt y ' 

A similar requirement will apply to the sixth order term considered here, providing another 
natural extension of the present work. 

In existing scenarios of Cold Electroweak Baryogenesis, the Higgs transition is triggered 
by a second scalar field. Although this would introduce model-dependence, including the 
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dynamics of this additional scalar field in the simulations is also an obvious next step (see 



also 



38 



40] for detailed simulations of such a system, but without CP- violation). 



For many years, Standard model CP-violation has been ruled out as a source for baryo- 
genesis. However, if baryogenesis was cold enough, so that the zero-temperature calculation 
of k cp can be trusted (or the similar calculation for CP-violating terms in the CP-even 
sector), it seems that the observed baryon asymmetry may in fact be accounted for without 
the addition of new sources of CP-violation. 
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Powers of D^cj) and An. 



Note that we use the <p upside down compared to the notation in 
Actually one for each factor, so (— l) 5 . 



Some of the literature (for instance 



30, 



31]) states the convention E® x = iTr[a a U n>x U^ x+0 ] 
while in fact using Eq. (|52p in practice. When used consistently, the two of course lead to 
equivalent results. 

A full-fledged determination of the k cp dependence is currently beyond our reach, numerically. 
The plaquettes are symmetrized forwards and backwards in space and time as for the action. 
Up to lattice errors of order 10 percent, and up to transients related to integer "flips". 
We checked that the same picture emerges when adding up in quadrature or absolute values 
of the force. 
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See for instance [37( for a simple 1+1 dimensional test. 
Which is absent in the Standard model 
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